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Abstract 



We estimate the masses of the 1 heavy four-quark and molecule states by combining exponential Laplace (LSR) and finite 
energy (FESR) sum rules known perturbatively to lowest order (LO) in a s but including non-perturbative terms up to the complete 
[dimension-six condensate contributions. We use double ratio of sum rules (DRSR) for determining the S U(3) breakings terms. The 
-S U(3) mass-splittings of about (50 - 1 10) MeV and the ones of about (250 - 300) MeV between the lowest ground states and their 
, „ 1 1st radial excitations are (almost) heavy-flavour independent. The mass predictions summarized in Table 2 are compared with the 
ones in the literature (when available) and with the three 7 C (4260, 4360, 4660) and 1^(10890) 1 experimental candidates. We 
conclude that the lowest observed state cannot be a pure 1 four-quark nor a pure molecule but may result from their mixings. We 
extend the above analyzes to the ++ four-quark and molecule states which are about (0.5-1) GeV heavier than the corresponding 
1 states, while the splittings between the ++ lowest ground state and the 1st radial excitation is about (300-500) MeV. We 
complete the analysis by estimating the decay constants of the 1 and ++ four-quark states. Our predictions can be tested using 
r ^_7 l ]some alternative non-perturbative approaches or/and at LHCi or some other hadron factories. 

Q^Keywords: QCD spectral sum rules, four-quark and molecule states, heavy quarkonia. 
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1. Introduction and a short review on the 1 ++ channel 
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A large amount of exotic hadrons which differ from the "stan- 
dard" cc charmonium and bb bottomium radial excitation states 
have been recently discovered in B-factories through _7i/ot + 7t~ 
and Yn + n~ processes and have stimulated different theoretical 
interpretations. Most of them have been assigned as four-quarks 
and/or molecule states [1]. In previous papers [2, 3], some 
of us have studied, using exponential QCD spectral sum rules 
(QSSR) [4-6] and the double ratio of sum rules (DRSR) [7- 
T0], the nature of the X(3872) 1 ++ states found by Belle [11] 
and confirmed by Babar [12], CDF [13] and DO [14]. If it is a 
(cq)(cq) four-quark or D - D* molecule state, one finds [2]: 

X c = (3925 + 127) MeV , with yff c = (4.15 + 0.03) GeV (1) 

corresponding to a ? c -value common solution of the exponential 
Laplace (LSR) and Finite Energy (FESR) sum rules. While in 
the _>-meson channel, using mt, - 4.26 GeV, one finds [2]: 

X b = (10144+ 104) MeV, with yfc = (10.4+0.02) GeV. (2) 

By assuming that the mass of the radial excitation X' Q * yftl, 
one can also deduce the mass-splitting: 



which is much lower than the ones of ordinary cc and bb states: 
i//(2S) - ift(lS) 590 ~ T(2S) - T(1S) 560 MeV, (4) 
suggesting a different dynamics for these exotic states. 

2. QCD Analysis of the 1" and ++ channels 

In the following, we extend the previous analysis to the case 
of the 1 and ++ channels and improve some existing analy- 
sis from QCD (spectral) sum rules in the 1 channel [15, 16]. 
The results will be compared with the experimental 1 can- 
didate states: 7(4260), 7(4360), 7(4660), 7,(10890). These 
states have been seen by Babar [17] and Belle [18, 19] and 
which decay into J/h/7t + 7t~ and Ytt + tt~ around the T(5S) mass. 
These states cannot be identified with standard cc charmonium 
and bb bottomium radial excitations and have been assigned to 
be four-quark or molecule states or some threshold effects. 

• Interpolating currents 

We assume that the 7 states are described either by the lowest 
dimension (without derivative terms) four-quark and molecule 

D S D* S vector currents given by: 



X' - X c a 225 MeV * X' - X h * 256 MeV 



(3) 
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• QCD input parameters 

The QCD parameters are given in Table 1 and we shall work 
with the running light quark parameters [20, 21]. 

Table 1: QCD input parameters. 



Parameters 



A(n f = 4) 

A(n f = 5) 

m s 

m c 

m b 

A, 

K = (SS)/{UU) 

M 2 

(a s G 2 ) 
3r$\ 



Values 



Ref. 



(324 ± 15) MeV 
(194 ± 10)MeV 
(0.114 + 0.021) GeV 
(1.26- 1.47) GeV 
(4.17 -4.70) GeV 
(263 ± 7) MeV 
(0.74 ± 0.06) 
(0.8 ± 0.2) GeV 2 
(7 ± 2) x 10~ 2 GeV 4 
(8.3±1.0)GeV 2 x<a s G 2 > 



p = (qqqq)l(qq) (2 ± l) 



[22-24] 

[22-24] 

[5, 24] 

[5, 24-27] 

[5, 24-26] 

[5] 

[9] 

[28-30] 

[22, 26,31-37] 

[26] 

[22, 28,31] 
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FESR 1.47 



' 5.0 5.5 6.0 6.5 7.0 7.5 8.0 
V»T [GeV] 



3. 1 four-quark state mass Yq 9 from QSSR 

In the following, we shall estimate the mass of the 1 four- 
quark state (Qq)(Qq), hereafter denoted by Yqa. In so doing, 
we shall use the ratios of the Laplace (exponential) sum rule 
and of FESR: 

^/(r)-M^-^f*. (7) 
• The Y cc i mass from LSR and FESR for the case b—0 
Using the QCD inputs in Table 1, we show the T-behaviour of 
M Ycd from 'R^d R m Fig- One can notice from Fig. la that 
the r-stability is obtained from yft[. > 5.1 GeV, while the t c - 
stability is reached for y/t^ = 7 GeV. The most conservative 
prediction from the LSR is obtained in this range of f c -values 
for m c = 1 .26 GeV and gives in units of GeV: 

4.79 < M Ycd < 5.73 for 5.02 < yJT c < 7 and m c = 1.26, 
5.29 < M Ycd < 6.11 for 5.5 < yfc < 7 andm f = 1.47. 

We compare in Fig. lb), the f e -behaviour of the LSR results 
obtained at the r-stability points with the ones from 'R F c ^ SR for 
m c =1.23 GeV (running) and 1.47 GeV (on-shell). One can de- 
duce the common solution in units of GeV: 

M Ycd = 4.814 for ^ c = 5.04(5) and m c = 1.26, 
= 5.409 for Vic = 5.6 and m c = 1.47. (8) 

We observe that the on-shell c-quark mass value tends to over- 
estimate Mji^ [3, 26]. The same feature happens for the eval- 
uation of the X{1 ++ ) four-quark state mass [2]. Then, we are 
tempted to take as a final result in this paper the prediction 
obtained by using the running mass m c (m c ) = 1262(17) MeV. 
Considering the uncertainties from the Table (1), we deduce 

M Ycd = 4814(57) MeV . (9) 

Using the fact that the 1st FESR moment gives a correlation be- 
tween the mass of the lowest ground state and the onset of con- 
tinuum threshold t c , we shall approximately identify its value 
with the one of the radial excitation. Assuming that, one can 
deduce the mass-splitting: M' Y ^ - M Y[d w 226 MeV, which is 
similar to the one obtained for the X(l ++ ) four-quark state [2]. 



Figure 1 : a) r-behaviour of My cd (1 ) from 'R^ R for the current mixing parameter b = 0, 
for different values of t c and for m c = 1 .26 GeV; b) The same as a) but for m c = 1 .47 GeV; 
c) / £ .-behaviour of the LSR results obtained at the r-stability points and comparison with 
the ones from H F cd SR for m c = 1 .26 and 1 .47 GeV. 

• The Ybd mass from LSR and FESR for the case b = 

We extend the previous analysis to the Z?-quark sector, in order 
to estimate the mass of four-quark state. Considering, like 
in the case of charm, as a final estimate the one from the running 
Z?-quark mass mt>(mb) = 4177(1 1) MeV [26], we deduce: 

M Ybd = 1 1256(49) MeV . (10) 

From the previous result, one can deduce the value of the mass- 
splitting between the 1st radial excitation and the lowest mass 
ground state: ML - M Yl , « M' v — M Y , ~ 250 MeV, which 

B r bd L U Y cd 'cd 

are (almost) heavy-flavour independent and also smaller than 
the one of the bottomium splitting (-560 MeV). 

In the following, we shall let the current mixing parameter 
b, as defined in Eq. (5), free and study its effect on the results 
obtained in Eqs. (9) and (10). In so doing, we fix the values of t 
around the r-stability point and t c around the intersection point 
of the LSR and FESR. The results of the analysis are shown in 
Fig. 2. We notice that the results are optimal at the value b — 0. 

For completing the analysis of the effect of b, we also study 
the decay constant f Ygd defined as: <0|^ 9 |y fid > = f YQd M YQd ^ ■ 
Doing the analysis, giving M Ygd and the corresponding t c ob- 
tained above, one can deduce the optimal values at b — 0: 

f Yci ~ 0.08 MeV and f Ybd 0.03 MeV , (11) 

which are much smaller than f„= 132 MeV, f p 215 MeV and 
/d — /b=203 MeV [38]. One can also note that the decay con- 
stant decreases like 1 /Mq which can be tested in HQET or/and 
lattice QCD. 

• 5 U(3) breaking for M Y<2s from DRSR 
We study the ratio M Yg JM YQd using DRSR: 

r i = where Q = c,b. (12) 



3 



a) 



-0.5 0.0 



b) 



-0.5 0.0 



Figure 2: a) /^-behaviour of My, for given values of r and t c and for m c = 1.26 GeV; b) 
the same as a) but for Afy, and for ny, = 4. 17 GeV. 

Extracting the values at r- and ^—stabilities points we can de- 
duce, respectively for \fT c = 5.1 and 1 1.6 GeV: 



r c sd = 1.018(1)^(5)^(2)^(2)^(1^, 

? sd = 1.007(0.5) mi ,(2) m ,(0.5),(l) 57 „(0.3) p . 



(13) 



Using the results for Yqj in Eqs. (9) and (10) and the values of 
the SU(3) breaking ratio in Eq. (13), we can deduce the mass 
of the Y Qs state in MeV: 



M Ycs = 4900(67) , 



M Ybs = 11334(55), 



(14) 



leading to the SU(3) mass-splitting: AM sd ~ 87 MeV 
AM 7 '] « 78 MeV, which is also heavy-flavour independent. 



The splitting (in units of MeV) with the first radial excitation 
approximately given by is: 

m j /4 ,s 2 ~ M JI*s 2 * 298 , M' TS2 - M TS2 * 213 . (19) 

In the same way, we analyse the t and t c behaviours of the 
S U(3) breaking ratios, from which, we can deduce: 

r* d = M ms JM ms , = 1.022(0.2) mc (5) ms (2)* , 

rj d = MrsJMr Sl = 1.01 l(l) m ,(2) m , (0.2), , (20) 

where S 3 = ss is a scalar meson. Then, we obtain in MeV: 



M 



= 5112(41) , M TS} = 10125(40) , 



(21) 



corresponding to the S U(3) mass-splittings: AM J J d ^ AM^ * 
1 10 MeV. Doing the same exercise for the octet current, we de- 
duce the results in Table 2 where the molecule associated to the 
octet current is 100 (resp. 250) MeV above the one of the sin- 
glet current for (resp T) contrary to the 1 ++ case discussed 
in [3]. The ratio of SU{3) breakings are respectively 1.022(5) 
and 1 .010(2) in the c and b channels which are comparable with 
the ones in Eq.(20). 



4. 1 molecule masses from QSSR 

• The D d( ^Dd( S ) and B^B^s) molecules 1 

Like in the previous case, we use LSR and FESR for studying 
the masses of the D d Dd and B d Bd and DRSR for studying the 
SU(3) breaking ratios: 

C = M D . D JM D:Dd , r* d = M B . A JM B * Bd . (15) 

Using the sets (m c = 1.26 GeV, -y/£ = 5.58 GeV) and (m b = 
4.17 GeV, V£= 11.64(3) GeV) common solutions of LSR and 
FESR, one can deduce in MeV: 

M D . Dd = 5268(24), Mg*g d = 11302(30) , 
rf d = 1.018(l) mt .(4)m s (0.8) ; ,(0.5k ( (0.2) p (0.1) C 3 , 
r B sd = 1.006(l) miJ (2)m. s (l),(0.5k ( (0.2) p (0.1) C 3 .(16) 

Using the previous results in Eq. (16), one obtains in MeV : 

M D . Ds = 5363(33) , M B ;b, = 1 1370(40) , (17) 

corresponding to a S U(3) mass- splitting: AM®f * 95 MeV * 
AM B J 68 MeV . 

• The J l>pS 2 and TS2 molecules 

Combining LSR and FESR, we consider the mass of the JI41S2 
and TS 2 molecules in a colour singlet combination, where S 2 = 
uu+dd is a scalar meson We work with the LO QCD expression 
obtained in [16]. Then we analyze the f c -behaviour of different 
r-extremas, from which we can deduce, for the running quark 
masses for V£ = 5.30(2) and 10.23(3) GeV, in MeV: 

M mS2 = 5002(31), M rSl = 10015(33) . (18) 



'Hereafter, for simplicity D and B denote the scalar D* and mesons. 



5. ++ four-quark and molecule masses from QSSR 

• Y°Q d mass and decay constant from LSR and FESR 

We do the analysis of the Y°, and Y?, masses using LSR and 
FESR match. We work with the current mixing parameter b — 
from which we deduce in MeV, for the running quark masses, 
and respectively for yfQ = 6.5 and 13.0 GeV: 

Myo = 6125(51) MeV , Myo = 12542(43) MeV . (22) 

cd bd 

One can notice that the splittings between the lowest ground 
state and the 1st radial excitation approximately given by yfQ 
is: M' - Myo « 375 MeV , M'„ - Mm =s 464 MeV , which 

is larger than the ones of the 1 states, comparable with the 
ones of the J/tfr and T, and are (almost) heavy-flavour indepen- 
dent. For completeness, we calculate the sum rule for the decay 
constants from which we deduce: 

fyo ^0.12 MeV and fa * 0.03 MeV , (23) 

cd bd 

which are comparable with the ones of the spin 1 case. 

• S f/(3) breaking for from DRSR 

Analyzing the t and t c behaviours of the S U(3), we deduce: 

r° s c d = 1.011(2) mi .(3.8)m. ! (1.4),(l) J7 „(0.7) p , 

r° s b d = 1.004(1)^(1.7^(0.3),, (24) 

leading to the masses in MeV and theirs respective S U{3) mass- 
splittings: 

Myo = 6192(59), Myo = 12592(50) , 

cs bs 

AM^, ^ 67 ^ AM 1 * 50 MeV . (25) 
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• M DdDd and M BdBd from LSR and FESR 

We consider as a final result the one corresponding to the run- 
ning masses, for ^T c = 6.25(3) and 12.02 GeV: 

M DdDd = 5955(48) MeV , M BdBd = 1 1750(40) MeV .(26) 

The splittings (in MeV) between the lowest ground state and 
the 1st radial excitation, approximately given by y7^, is: 

M' DdDd -M DdDd *290, M' BdBd -M BdBd *270, (27) 

which, like in the case of the 1 states are smaller than the ones 
of the J /if/ and T, and almost heavy-flavour independent. 

• S U(3) breaking for M-^ D ^ and M B B from DRSR 

Calculating the S U (3) mass ratios for the molecules D S D S and 
B S B S , considering different values of t c and the t c behaviour of 
their r-extremas, we deduce: 

r™ = M DsDs /M DdDd = lM5(l) mc (4)m s (2) K (l) liu (0.5) p , 
r °s! = M BsB JM BdBd = 1.008(l) mc (4)m. s (2),(l) I7 „(0.5) p . 

Using the previous values of M DdDd and M BdBd , we deduce: 
M DsDs = 6044(56) MeV , M BsBs = 1 1844(50) MeV, (28) 

which corresponds to a S U(3) splitting: 

AMff * 89 MeV * AMff * 94 MeV . (29) 



Table 2: Masses of the four-quark and molecule states from the present analysis 
combining Laplace (LSR) and Finite Energy (FESR). 



States 




States 




Four-quarks 


1— 




++ 


Y c d 


4818(27) 


cd 


6125(51) 


V 

1 cs 


4900(67) 


Y° 

1 CS 


6192(59) 




11256(49) 


1 bd 


12542(43) 


Ybs 


11334(55) 


bs 


12592(50) 


Molecules 


1~ 




++ 


D d D d 


5268(24) 


D d D d 


5955(48) 


D ; D S 


5363(33) 


D S D S 


6044(56) 


B d B d 


11302(30) 


BdB d 


11750(40) 


b" s b s 


11370(40) 


B S B S 


11844(50) 


Singlet current 


1~ 


Octet current 


1~ 


J/*s 2 


5002(31) 




5118(29) 


J/<fsS 3 


5112(41) 




5231(40) 


T5 2 


10015(33) 




10268(28) 


T5 3 


10125(40) 




10371(45) 



6. Summary and conclusions 



• The three F c (4260, 4360, 4660) 1 __ experimental candidates 
are too low for being pure four-quark or/and molecule DD* and 
J/1//S2 states but can result from their mixings. The F/,(10890) 
is lower than the predicted values of the four-quark and BB* 
molecule masses but heavier than the predicted YS 2 and YS 3 
molecule states. Our results may indicate that some other exotic 
structure of these states are not excluded. 

• For the 1~ , there is a regularity of about (250-300) MeV for 
the value of the mass-splittings between the lowest ground state 
and the 1st radial excitation roughly approximated by the value 



of the continuum threshold yf". at which the LSR and FESR 
match. These mass-splittings are (almost) flavour-independent 
and are much smaller than the ones of 500 MeV of ordinary 
charmonium and bottomium states. 

• There is also a regularity of about 50-90 MeV for the S U(3) 
mass-splittings of the different states which are also (almost) 
flavour-independent. 

• The spin states are much more heavier (> 400 MeV) than 
the spin 1 states, like in the case of hybrid states [5]. 
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